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ABSTRACT 

We study the mass-charge relation for the semiclassical extremal black hole 
of the S-wave sector Einstein- Maxwell theory coupled to conformal scalars. 
The classical ratio M/\Q\ = 1 is shown to be modified to M/\Q\ ~ 1 — k/Q 
for small k = Nh/ {12'kQ'^). Furthermore, numerical study for /c < 2 shows 
that M/\Q\ is a monotonically decreasing function of k. We speculate on the 
consequence of such a modification in the 4-dimensional context. 
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In connection with the physics of black hole evaporation, extremal black holes with 
vanishing temperature provide interesting theoretical laboratories. Immune to the Hawk- 
ing's thermal radiation, they are the first clues as to what the final stage of the evaporation 
process might be. But this does not mean we can consider the classical extremal black 
holes as the final product of the process. For one thing, the thermal behavior is already 
expected to break down for near extremal cases[]T]]. Zero Hawking temperature simply 
means the leading quantum effect disappears. In order to address the questions of black 
hole quantum physics, we need a more systematic way of treating quantized matter in 
nontrivial geometries. In two-dimensional models, such a method has been adopted by 
Callan et al. (CGHS)0, and used extensively to study 2-D black holes semiclassically. 

In this letter, we want to concentrate on the case of the extremal Reissner- Nordstrom 
black hole and to study how semiclassical effects modify one of the classical properties, 
namely ADM mass M. The model we consider is dimensionally reduced Einstein-Maxwell 
theory. By restricting to the spherically symmetric sector we obtain the following 2-D 
act ion ,0 

Sg = -^j dPx v/^e-^-^ + 2 (V0)' + 2e^^ - F^), (1) 
where the 4-D metric is split into 2-D metric g^"^^ and the dilaton part 

g(^)=gi^)+e-'Un'. (2) 

The finite mass solutions with regular horizons are the well-known Reissner-Nordstrom 
solutions with mass M and charge Q satisfying the inequality M > \Q\. 

g(^) = _F(r) dt' + ^ + dn\ F(r) = 1 - ^ + ^. (3) 

When the inequality is saturated, F{r) has a double zero at the horizon r = M = 

\Q\ and the corresponding extremal black hole has zero Hawking temperature, hence no 

thermal radiation emanates from the horizon a long time after the black hole formation. 
■^G = c = 1 in this letter 
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This implies that the usual late time estimate of the Bogoliubov transformation is not 
the leading quantum correction. It vanishes identically and we need to study next the 
nonvanishing contribution, which may or may not depend on the history of the collapse. 

For this purpose, we can follow CGHS and couple conformal scalars to the above 
2-D action. One can regard these 2-D scalars as the S-wave part of massless 4-D fermions, 
alternatively. Integrating them out completely, which is possible since we are in a two- 
dimensional toy world, produces the non-local Polyakov-Liouville action with a partic- 
ular coefficient, which summarizes the effect of the quantized matter on gravity. Further- 
more, this effective semi-classical gravity can be conveniently handled with the introduc- 
tion of a scalar field z with a background charge in the following manner Q], 

S = S,-^Jd'x ^ {{Vzf - zR^^^). (4) 

Since the field equation for z reduces the second term to the original Polyakov-Liouville 
action of central charge A^, solving this theory at tree level is equivalent to studying the 
semi-classical theory of 2-D gravity coupled to scalars. These are the leading terms 
in the large N expansion of the full 2-D quantum theory, where is large but Nh is of 
order one. 

Notice that the scalar curvature i?*^^-* acts as an external source coupled to the z field. 
This effectively induces the usual Hawking radiation in a classical black hole geometry, 
i.e., a stationary point of Sg only. Because the classical black holes radiate, the theory does 
not have any static solution with finite mass and regular horizon of finite temperature. 
The only static solutions of finite mass are those of zero temperature, which were first 
studied by S. Trivedi0. 

But first let us consider the effect of quantizing the original A^ conformal scalars in 
a classical background. For example, given a classical geometry, the expectation value 
of the matter energy momentum tensor can be found simply by evaluating the classical 
energy momentum tensor of 2;-field on that classical background. A family of classical 
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geometries known as the Vaidya metric is particulary relevant to our discussion. 

g^'^ = -(1 - + ^)dv' + 2dvdr + d^f (5) 

It represents a collapsing massless shell whose cumulative energy and charge at retarded 
time V are m{v) and e{v). For smooth m and e^, the cosmic censorship is achieved by 
requiring the positive energy condition for the shelipj.0 
For our purposes, however, it is appropriate to choose 

m{y) = Me{v -Vq), e^{v) = Q'^Oiy - Vq), (6) 

where 9 is the usual step function. The geometry is then that of an initial Minkowski 
spacetime glued to a Reissner-Nordstrom black hole across an ingoing null shock wave 
located at f = vq. Introducing a new coordinate u = v — 2 J F~^{r) dr with F{r) as in 
eq. (^, {v,u) form a pair of light-cone coordinates above the shock, 

g{2) ^ _p(^r) dvdu, V > vq. (7) 

In these coordinates, f — oo is the future null infinity, and m — ^ oo is the future event 
horizon. Suppose we impose an initial condition on the N matter fields such that the 
expectation value of the energy-momentum tensor vanishes in the Minkowskian region. 
Using energy-momentum tensor conservation, this can be translated into 

Nh 

< Tun >\v=vo = {-TTridlp - (dupf) + tuuiu)) \v=vo= (8) 

P = (^logi^), 

where tuu comes from the homogeneous part of the solution to z field equation. This 
implies the following form of < Tuu > as we approach the future null infinity. 

< Tuu >|™= -tuuiu) = ^(^^'W - lF{r)F"{r)) \u=uo-2 J f-^ dr ■ (9) 
^Of course, we have introduced an external charged matter source to create the shell itself. 
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As M oo, this clearly shows a steady flux proportional to the temperature squared 
(F'(r — > horizon) ~ Tuawking)- Also as expected, this asymptotically steady flux is 
absent, if M is equal to \Q\ so that F'{r horizon) = 0. However, there is a finite 
integrated flux; the total energy radiated is 

AM= / <Tuu>\v^oodu= — dr = -\Ql (10) 

J-oo 9o7r J\Q\ r 

where we used F' {r = \Q\) = Q for the M = \Q\ case. Since we ignored the gravitational 
backreaction, this estimate is valid only for small k = Nh/{12TTQ'^), or equivalently for 
large black holes. Notice that AM is positive for any F > with a double zero at the 
event horizon. AM represents the energy radiated away by the quantized matter, and 
after properly taking into account the gravitational backreaction, the Bondi Mass of the 
system should approach as u — > cxo 

|g|(i-^ + o(fc2)). (11) 

One might assert that it is not clear whether the estimated loss depends on the particular 
history of the collapse chosen. After all, the metric chosen can never be realized, since one 
cannot assimilate the collapsing process by a smooth version of the shock wave. As shown 
in 1^, for smooth m{v) and e^(f) satisfying the positive energy condition, the extremality 
can never be achieved in finite time. It might be that a realistic collapse scenario produces 
different AM. We will show that the above estimate of energy loss is robust by finding 
numerically the ADM mass of the semi-classical analogue of the extremal black hole which 
must be the end stage of the processs described so far. 

Semi-classical static solutions with extremal horizons have been studied near the 
horizon^]. The requirement of zero temperature horizon specifies a unique initial condi- 
tion at the horizon for a given total charge, and the resulting static solution is known to 
be asymptotically flat. There is no known analytical form of the solution, but it is, in 
principle, possible to carry out numerical integration. 

Before going into details of the simulations performed, it is helpful to discuss other 
static solutions of finite mass, all of which have naked singularities. Those with smaller 
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masses, to be called supercritical, are qualitatively similar to the classical ones with 
M < Q. The radius e""^ monotonically decreases as we approach the naked singularity 
at near origin. On the other hand, solutions with larger masses, to be called subcritical, 
are quite different from classical analogues M > Q, which have curvature singularity 
at the center e"'^ — hidden by two layers of nonextremal horizons, since we assume 
no heat bath to support nonextremal horizons. (Heat bath makes ADM mass infinite.) 
More specifically, a semi-classical subcritical solution has a lower bound on the value of 
the radius e~'^ near would-be horizon. One can distinguish the two species by observing 
whether the simulation stops in the middle or continues all the way to the critical value 
of the radius e~'^'='' = ^JkQ'^. 

Coming back to the actual simulation, it turns out that static field equations can be 
decoupled to produce a single first order differential equation with the following gauge 
choice. 

In this gauge we can extract two independent first order differential equations. 



*4)^ + 2r4) + (l-B'^+^) = 

Solving for A' /A in terms of B produces a first order differential equation for l/B"^. We 
performed two independent simulations. First, we started from the asymptotic region 
with the initial condition determined by M/\Q\, and searched for the range of M/\Q\ 
producing an extremal horizon (a double zero of 1/B'^). Secondly, we integrate outward 
from the known behavior near the extremal horizon and extract the ADM mass by fitting 
the curve in the asymptotic region. Since the initial points are near, but not quite at the 
horizon or r = oo, we needed to calculate accurate initial conditions. Symbolic expansions 
of 1/B^ in appropriate coordinates, solving the equation above approximately, are used 
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for this purpose. Fortunately, the nonanalytic behavior of the metric near the horizon 
emphasized in does not occur for 1/5^ as a function of r. We used MATHEMATICA 
for all numerical and symbolic calculations as well as preparation of the plot. As we 
improved the accuracy of the numerical calculation by supplying more accurate initial 
data, and also by increasing the intrinsic accuracy of the program used, the results from 
each simulation converge to each other. The data for M/\Q\ obtained by the two methods 
coincide to an accuracy of 10~^. 

The simulation is carried out only for k < 2 because the extremal horizon disappears 
beyond k = 2, when the horizon radius is equal to the critical value of the dilaton 
e-<^- = ^/kQ^. The plot of M/\Q\ as a function of k = Nh/{127iQ^) (Fig. 1) clearly 
shows the initial slope of —1/6 calculated above. Furthermore, up to /c = 2, the ratio 
continues to drop as we increase Nh or decrease the charge \Q\. 

What can we learn from this little demonstration? The first and foremost fact is that 
higher order corrections to Hawking's calculation must be taken account into even for such 
a crude operation as mass measurement. One should expect that a similar mechanism 
works for four-dimensionally black holes and the classical bound M > \Q\ is modified, 
unless some unbroken extended supersymmetry protects it. But the model we used gives 
few clues as to what the modification might be. While 2-D conformal scalars can be 
interpreted as the S'-wave modes of 4-D massless fermions, we cannot regard our model 
as a quantitative approximation to the full 4-D physics. There is no generic mass gap 
present to separate S'-wave fermions out from the rest. 

Nevertheless this doesn't prevent us from speculating on the effect of such a modified 
mass-charge relation in 4-D. In particular, suppose the same monotonic decreasing behav- 
ior is realized for the four-dimensionally extremal black holes. The possibility has been 
contemplated by J. Preskill for electrically charged extremal black holes with emphasis 



on charge renormalization ||10| . The most immediate consequence would be to lift the 
well-known degeneracy for multi-extremal black hole configurations. Classically, a family 



7 



of solutions known as the Papapetrou-Majumdar space-time describes many extremal 
black holes at rest relative to one another. The total ADM mass of such a solution is the 
sum of the individual masses, 

M = Y.m- (14) 

i 

This can easily be seen by imagining each hole separated from one another far away, so 
that whatever potential energy there might be becomes negligible. In fact, there is no 
potential between individual black holes and the total mass is given by eq. (|T^) for any 
finite separations. Therefore two different multi-black hole configurations in equilibrium 
have the same energy provided that the sum of absolute value of the charges are equal. 
But, with the modified M/\Q\ which decreases as decreases, the same reasoning shows 
that it is energetically favorable to split one big black hole into many smaller ones. The 
classical degeneracy is lifted. 

Classical physics forbids such a bifurcation process, since it violates the second law of 
black hole thermodynamics. However, there has been suggestions of possible finite action 
instantons mediating bifurcation of the extremal Reissner-Nordstrom black holes. In fact, 
D. Brill found an instanton of finite action interpolating between two Bertotti- Robinson 
metrics with different numbers of necks [IT]. It is well known that a Bertotti- Robinson 
metric with a single neck approximates an extremal Reissner-Nordstrom black hole near 
the horizon. If the initial and the final states are of the same energy, the instanton will 
take infinite Euclidean time to make the transition, and the stationary state would be a 
linear combination of the two classical configurations. With the modified M/\Q\ relation 
however, a relevant Euclidean solution is a bounce solution and a big extremal Reissner- 
Nordstrom black would decay to many extremal black holes of smaller charges distantly 
separated.0 

So far, we have completely ignored the possible presence of charged matter fields. 
Suppose there is an elementary charged particle of mass m and charge e and consider 



similar observation has been made in the context of classical dilatonic black holes with a massive 
dilatonill. 
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an extremal black hole of mass M and charge Q. For m « |e|, the Schwinger pair 
production near the horizon is always dominant over a possible bifurcation process and the 
black hole charge will eventually be wiped out. But for sufficiently large m > rrimin, it will 
be kinematically impossible for an extremal black hole to lose its charge by emiting these 
charged particles [|10|] . For a large black hole M ^ m in particular, we have mmj„/|e| ~ 
M/\Q\. Therefore the model we considered should be regarded as a possible scenario for 
magnetically charged extremal black holes in a world where the magnetic monopole comes 
with mass comparable to, or even larger than, its charge in Planck units. 

In summary we have found that the classical inequality M > \Q\ can be modified 
through semi-classical effects. It would be most interesting to find out about similar 
effects in the context of 4-dimensional models but it is beyond the scope of this letter. 

We would like to thank J. Preskill for critical reading of the manuscript. Also P.Y. 
thanks S.Trivedi for invaluable discussions at various stages of this work. 
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Figure Caption 

Fig.l: Plot of M/\Q\ versus k = Nh/{127iQ^). The straight hne shows the leading 
behaviour M/\Q\ = 1 — k/6. The dots arc the actual numerical results from the two 
independent simulations. Data points are at k = n/10 for n = 1, 19 as well ask — 0.001. 
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